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Dynamical crossover in water is studied by means of computer simulation. The crossover temper¬ 
ature is calculated from the behavior of velocity autocorrelation functions. The results are compared 
with experimental data. It is shown that the qualitative behavior of the dynamical crossover line 
is similar to the melting curve behavior. Importantly, the crossover line belongs to experimentally 
achievable (P, T) region which stimulates the experimental investigation in this field. 

PACS numbers: 61.20.Gy, 61.20.Ne, 64.60.Kw 


In recent decades wide attention of researchers is at¬ 
tracted to the field of supercritical fluids. Supercritical 
fluids are of extreme importance both for fundamental 
research and for practical applications especially in chem¬ 
ical industry. In this respect it is important to give an 
unambiguous definition of the supercritical state itself. 
For some years a so called Widom line, which is a line of 
supercritical maxima of correlation length and thermo¬ 
dynamic response functions (isobaric heat capacity cp, 
isothermal compressibility fcp, heat expansion coefficion 
ap) in fluids, was used to extend a liquid-gas coexistence 
line into supercritical region [i|. Later on it was shown 
that Widom line is ill defined and can not be used for de¬ 
marcation of a supercritical region of fluid ii. Several 
authors reported Widom lines of different systems and 
arrived to the same conclusion . One can clearly see 
that another way of demarcation of gaslike and liqnid-like 
fluids beyond the critical point was necessary. 


Such a way was proposed in our previous publications 
iSil- These publications introduce a so called Frenkel 
line which is the line of dynamical crossover in fluids. 
It was shown that below the Frenkel line the particles 
of fluid make few oscillations at some quasi-equilibrium 
position following a jump to another quasi-equilibrium 
point. This model was proposed by J. Frenkel 12| after 
whom the line was named. Above the Frenkel line the 
particles of fluid move like in a gas, by long jumps before 
a collision with another particle occurs. As a result be¬ 
low the Frenkel line liquid properties demonstrate some 
solidlike behavior whilst above it the properties of liq¬ 
uid are similar to a dense gas ones. Impact of crossover 
from crystallike microscopic dynamics of fluid particles 
to a gaslike one on different properties of the fluid was 
addressed in details in Refs. @,[13 • In Ref. 11| it was 
shown that the most convenient way to find the location 
of Frenkel line in the phase diagram is by monitoring a 
velocity autocorrelation function (vacf) of the fluid. Bas¬ 
ing on this criterion Frenkel line of several model systems 
(Lennard-Jones and soft spheres Iffl) and realistic ones 


(liquid iron [iSLjCarbon dioxide [1, [I^, r/P4P/2005 
model of water 1^, methane [l^ and hydrogen [l^) was 
calculated. This work extends the investigations in the 
field to the most important liquid - water. 

Phase diagram of water is extremely complex. It con¬ 
tains numerous solid phases including both crystalline 
and amorphous ones. However, the fluid part of the phase 
diagram is no less interesting. In addition to usual liquid 
- vapor transition it is widely assumed that a liquid-liquid 
phase transition (LLPT) takes place in unachievable re¬ 
gion of (P,T) parameters [l^. Water demonstrates a set 
of liquid state anomalies such as density anomaly, diffu¬ 
sion anomaly, structural anomaly and many others 171. 
Moreover, water has multiple Widom lines 18[. The first 
Widom line is related to liqud-gas transition while the 
second one can be attributed to the hypothetic LLPT. 
Both Widom lines of water were vividly discussed in liter¬ 
ature (see, for example JM, 3“ 231 and references therein). 
In particular, in Ref. [22| a connection between Widom 
line and dynamical properties of water was proposed. 
Since Frenkel line is the line of dynamical crossover in 
fluids some kind of relation between the Widom line and 
Frenkel line can exist. 

In the present work we study the behavior of wa¬ 
ter by means of molecular dynamics simulations. An 
SPC/E model of water is used [ 23 . The phase dia¬ 
gram of this model was reported in several publications. 
In a comparison of solid part of the phase diagram 
of SPC/E, several variants of TIPAP model and ex¬ 
perimental results is given. One can see that all mod¬ 
els fail to reproduce the whole complexity of the exper¬ 
imental phase diagram, but manage to describe some 
parts of it. In particular, SPC/E model is good in 
describing boiling curve of water. In Ref. boiling 
curve of SPC/E water is reported. The critical param¬ 
eters are found to be Tc = 651.TAT, pc = 0.326g/cm^ 
and Pc = 189bar. Experimental critical point of wa¬ 
ter corresponds to Tc = GdT.lSiF, pc = D.iTlg/crv? and 
Pc = 220.556ar. One can see that except some difference 
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in critical pressure the critical point of SPC/E model is 
very close to the experimental one. Moreover, in Ref. 


27[ a comparison of SPC/E model and ab-initio re¬ 


sults at high pressures and high temperatures was re¬ 
ported. It was shown that the discrepancy of SPC/E 
model and ab-initio results is of the order of 15 — 20% for 
T = lOOORT and pressure up to about IQOkhar. However, 
at T = 2000itr and pressures up to approximately 90kbar 
the agreement of ab-initio and SPC/E results is within 
5% which should be considered as good agreement. One 
can guess that at T = lOOORT and so high pressure the 
results are affected by crystallization effects while at tem¬ 
peratures well above the melting line SPC/E model can 
be used to study the high pressure behavior of water. 

A system of 4000 water molecules in a cubic box 
was simulated in molecular dynamics at constant vol¬ 
ume, number of particles and temperature (canonical en¬ 
semble). The temperature was held constant by Nose- 
Hoover thermostat. The density was varied from pmin = 
Q.ig/cwS' up to Pmax = 2.0g/cm^ and the temperatures 
from Train = 275K up to Tmax = 5 ■ 10^AT. Initially 
the system was equilibrated for Ins with a time step 
dt = 1/s. After that it was simulated more Ips with the 
same time step in order to calculate the thermodynamic 
properties. Finally, 10® steps with timestep dt = 0.1/s 
were made in order to well reproduce the decay of veloc¬ 
ity autocorrelation function (vacf). 

All simulations were performed using lammps simula¬ 
tion package (28| . 

As it was proposed in our earlier publications several 
methods to find the location of Frenkel line can be used 
The most convenient one is based on the lose of 
oscillations of vacfs. This criterion is used in the present 
work. Figs. [T] (a) and (b) show the vacfs of oxygens for 
two densities: p = 1.0 and 1.3 g/crn?. One can see that 
the low temperature vacfs for these two densities look 
qualitatively different while at high temperatures they 
become very similar. In particular, as the temperature 
increases the oscillations of vacfs become less pronounced 
and finally disappear. 

One more way to estimate the location of Frenkel line 
in P — T or p — T diagram is related to heat capacity 
of liquid |9l4llj. In case of water the heat capacity cy 
undergoes strong decay upon isochoric heating. Next to 
the melting line the heat capacity per molecule is about 
9kB or Sfcs pet atom is Boltzmann constant) while at 
high temperatures it becomes as small as l.Sfcs per atom. 
The heat capacity per molecule at Frenkel line should 
be 6kB per molecule or 2kB per atom. The location of 
Frenkel line by cy criterion was evaluated from experi¬ 
mental data. The data were taken from NIST database 


37|. 

Frenkel line of water is shown in Figs. [2] (a) and (b). 
From Figs. [2] (a) one can see that Frenkel line starts at the 
boiling curve at temperature Tp = dbOAT « 0.7Tc, where 
Tc is the critical temperature. In our previous publica- 
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FIG. 1: Normalized velocity autocorrelation functions of oxy¬ 
gen atoms at density (a)p = l.Og/cm^ and (b) p — 1.3g/cm^. 
The inset in panel (b) enlarges the time period 40 — 140 ps 
where the oscillation behavior of vacfs takes place. (Color 
online). 


tions it was shown that the same ratio Tp/Tc takes place 
in Lennard-Jones fluid and liquid iron. 

One observes extremely fast grow of the Frenkel line 
temperature in a narrow interval of densities p = (1.2 — 
1.22)g/cw?. It is related to extremely slow disappearance 
of vacf oscillations in this region. 

Figs. [H (b) shows the location of Frenkel line of wa¬ 
ter in P — T phase diagram. Phase diagram of water is 
very complex. It demonstrates numerous solid phases. In 
particular, usual ice I melting line has a negative slope. 
Other solid phases melting lines have positive slopes. The 
shape of the Frenkel line of water qualitatively resembles 
the shape of the melting line. At low pressures Frenkel 
line very slowly increases with pressure. Later on at pres¬ 
sure about lOkbar the slope of the Frenkel line like the 
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FIG. 2: Location of Frenkel line of water in the phase diagram 
in (a) p — T and (b) P — T coordinates. In case of P — T 
diagram the boiling curve corresponds to the experimental 
curve [ 2 ^, the critical point is taken for SPC/E model [^. 
The melting line is combined from several publications [30l - 
1^ . Experimental data for heat capacity are taken from NIST 
database f^. The Frenkel line of T7P4P/2005 water is taken 
from Ref. Il'^ - Widom line of TIP4P/2005 model is taken 
from Ref. [^. Panel (c) enlarges the moderate pressure part 
of the diagram. The notation of this panel is the same that 
in panel (b). (Color online). 


one of the melting curve rapidly increases. 

Importantly, the relation between the temperature at 
the Frenkel line Tp and the melting temperature Tm 
changes upon increasing the pressure. At low pressures 
(before the rapid increase) the ratio Tp/T^ is close to 
2. At pressures P lOOkbar this ratio increases up to 
5. On further rise of pressure it reaches the value of 9 
at P 900fc6ar. It means that in the range of pressures 
considered in the present work the Frenkel line bends up 
with respect to the melting line. 

In our previous publications it was proposed that in 
the limit of high pressures Frenkel line should be parallel 
to the melting line in double logarithmic coordinates. In 
case of water we are not aware of any measurements of 
the melting curve above 1000A:6ar. We expect that the 
Frenkel line and the melting curve will be parallel in the 
high pressure limit, but one needs to extend the melting 
curve to higher pressures in order to check it. Although 
this conclusion may be violated by transition of water 
into superionic phase under high pressure . This phe¬ 
nomena can not be taken into account in frames of purely 
classical model used in the present work. 

Very recently Frenkel line of water calculated by vacf 
criterion for a different model (T/P4P/2005) was re¬ 
ported 0. This line is shown in Fig. [2] (b) for the 
sake of comparison. One can see that this line is sys¬ 
tematically higher then our line. However, the lines are 
very close to each other and this small discrepancy can 
be attributed to the different models under investigation. 
The authors of 0 studied the Frenkel line of water up to 
Pmax = SOkbar. In our work the Frenkel line is traced up 
to pressures as high as almost Pmax = 2000kbar where al¬ 
lowed us to see an interesting phenomenon. From Figs. [2] 
(b) one can see that the bend of Frenkel line in dou¬ 
ble logarithmic coordinates takes place at pressure about 
50kbar. One can relate this bend to some changes in the 
local structure of the liquid. 

Figs. [2] (b) and (c) show also the location of Widom 
line of T/P4P/2005 model of water calculated from max¬ 
ima of isobaric heat capacity cp in Ref. [ 2 ^. As it 
was shown in several recent publications the supercritical 
maxima of different substances rapidly vanish on depart¬ 
ing from the critical point In Ref. the Widom 

line extends up to the pressure O.SSkbar which is lower 
then the Frenkel line starts. Moreover, extrapolation of 
the Widom line to higher pressures should go above the 
Frenkel line of both SPC/E and T/P4P/2005 models of 
water. From this one can conclude that these lines are 
not related to each other. One can assume that at low 
pressure the crossover of dynamical properties of liquid 
is governed by the Widom line while at higher pressures 
it is determined by the Frenkel line. 

Let us consider the coordination number of oxygens. 
The coordination number can be calculated from radial 
distribution function (rdf): NN = iirp g{r)r^dr, 
where p is the number density and Cmm is the location of 
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FIG. 3: Number of nearest neighbors of oxygen along T = 
lOOOif isotherm. For comparison literature data of NN along 
the melting line are also given. The data are taken from Ref. 
[ 3 ^ (Week) and Ref. [s^ (Katayama). (Color online). 


the first minimum of rdf. The number of nearest neigh¬ 
bors (NN) along the isotherm T = lOOOitT is shown in 
Fig. [3l NN of water along the melting line was reported 
in Refs, [sl, [s^. Although our results belong to an 
isotherm while the results of these publications are re¬ 
lated to the melting line, the NNs are calculated in sim¬ 
ilar pressure interval. That is why we show them in the 
same plot for comparison. The difference between this 
work and the literature data should be referred not only 
to different lines in P — T plane but also to different mod¬ 
els studied and different methods of calculation. One 
can see that up to pressure of approximately bOkbar the 
coordination number rapidly increases while above this 
threshold it holds approximately constant. One can con¬ 
clude, that at small pressures the local structure is very 
sensitive to the pressure change while at higher ones the 
local structure is very stable which is similar to the case 
of simple liquid. Therefore one can say that water be¬ 
comes ’’simpler” upon increasing pressure. 

We expect that similar rapid increase of Frenkel line in 
the density range p = 1.2 — 1.22g/cm^ can be observed 
in T/P4P/2005 as well. However, the results reported 
in the work [l^ end up at the density 1.2g/cm^. So one 
needs to extend them in order to check this assumption. 

We report a computational study of dynamical 
crossover in water. The temperature of crossover (the 
Frenkel line temperature Tp) is calculated from the be¬ 
havior of velocity autocorrelation functions. The results 
are compared to the experimental ones obtained from 
isochoric heat capacities. It is shown that qualitative be¬ 
havior of Frenkel line is similar to the behavior of the 
melting curve. However, the ratio Tp/Tm increases with 
increasing of pressure. It means that the lines diverge 


in the range of pressures considered in the present work. 
This divergence can be related to the change of the local 
structure of water upon increasing the pressure. 

Importantly, the temperatures of dynamical crossover 
appear to be rather low at moderate pressures. At pres¬ 
sures P < 30 kbar Tp does not exceed 1000 K. It means 
that {P,T) parameters of Frenkel line can be achieved 
experimentally. Such experimental works would be im¬ 
portant not only for deeper understanding of dynamical 
behavior of liquids but also could serve for supercritical 
technology. Importantly, the properties of fluids at the 
Frenkel line are close to the optimal ones for technologi¬ 
cal applications [i^ . Water as well as carbon dioxide are 
amon g th e most important and widely used supercritical 
fluids [4lj | . Supercritical water is used for many different 
applications, such as green solvent, as reaction medium 
for different chemical processes, for production of biofuel, 
oxidation of hazardous materials which is important for 
dangerous waste disposal. Applications of supercritical 
water include separation, extraction and purification of 
different substances and many others ^|. 

The principal property of supercritical water providing 
its widespread application is its solving power. In Ref. 
[ 1 ^ the solubility of different solutes in carbon dioxide 
and its relation to the Frenkel line was discussed. It was 
shown that the solubility maxima are close to the Frenkel 
line. However, the data for solubility maxima of differ¬ 
ent solute in supercritical water are not available at the 
moment. Basing on the discussion of 1^ and the re¬ 
sults of the present paper one can propose that the opti¬ 
mal solving power of water should belong to the interval 
T = 700 - IOOOA: and P = 1 - 3GPa. Currently, the 
most widespread usage of supercritical water belongs to 
the temperatures interval 650 — lOOOAT and to the pres¬ 
sures up to O.bkbar. However, these {P,T) conditions 
were found empirically and do not have any solid the¬ 
oretical ground. Moreover, up to now the supercritical 
technology advances the theoretical foundations in the 
field. This work as well as Ref. 1^ allow to predict the 
best (P, T) conditions for supercritical water application 
which makes these publications the pioneering works in 
developing the theoretical basis of supercritical technolo¬ 
gies. 
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